###5.1 Editando alguns parmetros grficos
##5.Crie um grfico de disperso entre dap1) e ht2) com: 

# Legendas dos eixos com nomes das variveis e suas unidades
# Marcaes do eixos (ticks) para dentro da rea do grfico 
# Apenas dois eixos (formato L)
# Ttulo informativo
# Tamanho das fontes maiores que o padro

e.saligna<-read.csv("esaligna.csv")
str(e.saligna)
head(e.saligna)
summary(e.saligna)
names(e.saligna)

#nomes
dap<-"Dimetro  altura do peito(cm)"
tit.saligna<-"Relao DAP x Altura"

par(cex=1.2)
plot(e.saligna[,4]~e.saligna[,5], xlab="Altura(m)",ylab=dap,
          tcl=0.4,bty="l",main=tit.saligna)

##########################
###5.2 Dois grficos juntos

# Use as variveis dap e talhao para construir dois grficos, 
#colocando-os lado a lado. O primeiro deve ser um grfico de desenho de 
#caixa (boxplot) da varivel dap em funo do fator talho. 
#O segundo deve ter apenas a mdia e uma barra de desvio-padro do dap, 
#para cada talho.
# Insira tambm uma letra para dizer qual  o grfico a e qual  o b
# (tanto faz, quem  um e quem  outro).

par(mfrow=c(1,2))
names(e.saligna)
as.factor(e.saligna[,3])

#primeiro
boxplot(e.saligna[,4]~e.saligna[,3],ylab="DAP(cm)",xlab="Talho")
text(6.5,23,"a")

#segundo
media.dap<-tapply(e.saligna[,4],INDEX=e.saligna[,3],FUN=mean)
desvio.dap<-tapply(e.saligna[,4],INDEX=e.saligna[,3],FUN=sd)

(dap.talhao<-matrix(c(media.dap,desvio.dap),nrow=6,ncol=2))
rownames(dap.talhao)<-c(16,17,18,22,23,32)

desv.pos<-dap.talhao[,1]+dap.talhao[,2]
desv.neg<-dap.talhao[,1]-dap.talhao[,2]

plot(dap.talhao[,1],xaxt="n",ylab="DAP(cm)",ylim=c(0,22),xlab="Talho")

axis(side=1,labels=c(16,17,18,22,23,32), at=1:6)
text(6,22,"b")

arrows(seq(1:6),dap.talhao[,1],seq(1:6),desv.pos,angle=90, code=2,length=0.05)
arrows(seq(1:6),dap.talhao[,1],seq(1:6),desv.neg,angle=90,code=2,length=0.05)

#####################

#### 5.3Adivinhando o cdigo
##Figura 1

ex.3<-read.csv("exercicio3.csv")
str(ex.3)
head(ex.3)
summary(ex.3)
names(ex.3)

par(mfrow=c(1,2),mar=c(5,4,2,1))
ED<-"Euclidean distances"

par(bty="L",tcl=0.3,pch=17,cex=1,cex.axis=0.8)
plot(ex.3[,1],ex.3[,2],xlab="Log(Patch size)(ha)",ylab=ED,ylim=c(0,3.0),
       xlim=c(0.5,2.0),cex.lab=1.1)

#model.ex.3 <-  lm(ex.3[,1]~ex.3[,2])
#abline(model.ex.3)#no  identica# intersepto e inclinao diferentes

abline(2.5,-0.68,lwd=2.5)

text(2.05,3,"a")

###Figura 2

ex.3$y2<-as.factor(ex.3[,4])
class(ex.3[,4])

plot(ex.3[,4],ex.3[,3])
identify(ex.3[,4],ex.3[,3])

#27- linha que esta o ponto,excluir essa linha
#no deu certo excluir linha pq boxplot muda a forma

plot(ex.3[,4],ex.3[,3],ylim=c(0,3.0),,xaxt="n",tcl=0.3)

plot.box<-c("","Medium\nEdge","Medium\nInterior",
            "Large\nEdge","Large\nInterior","")

mtext("Small",side=1,at=1,line=0.5,cex=0.9)
axis(side=1, at=1:6,label=plot.box)
mtext("Control",side=1,at=6,line=0.5,cex=0.9)

#inserir figura pra ficar acima do ponto circular (ponto da linha 27)
rect(4.7,1.6,5.2,1.8,col="white",border=NA)

text(1:5,3,c("*","*","**","*","***"))
text(6.5,3,"b")



##################         ##############################